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Tapering studies for Terawatt level X-ray FELs with a superconducting undulator 
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We study the tapering optimization scheme for a short period, less than two cm, superconducting 
undulator, and show that it can generate 4 keV X-ray pulses with peak power in excess of 1 terawatt, 
using LCDS electron beam parameters. We study the effect of undulator module length relative to 
the FEL gain length for continous and step-wise taper profiles. For the optimal section length of 
1.5m we study the evolution of the FEL process for two different superconducting technologies NbTi 
and NbSSn. We discuss the major factors limiting the maximum output power, particle detrapping 
around the saturation location and time dependent detrapping due to generation and amplification 
of sideband modes. 

PACS numbers: 41.60.Cr, 41.60.Ap 


I. INTRODUCTION 

Much recent scientific effort has been devoted to study¬ 
ing the possibility of using a tapered undulator [1] to 
achieve terawatt level hard x-ray pulses in the next gener¬ 
ation of Free Electron Lasers (FELs) [2] [3] [1] [5] [S] . The 
motivation for pursuing this goal comes primarily from 
the bioimaging community where a terawatt power coher¬ 
ent x-ray source would open the doors to single molecule 
imaging mEM- Recent numerical work m E] shows 
that a self seeded hard X-ray FEL with LCLS-II like pa¬ 
rameters in a 200-m permanent magnet undulator has the 
capability of reaching TW level pulses with the longitudi¬ 
nal and transverse coherence necessary for coherent X-ray 
imaging [12) . Spatial constraints, an increase in tunabil- 
ity and a longer machine lifetime have since led towards 
considering using superconducting technology for the un¬ 
dulator design rather than permanent magnets Ba¬ 
in this study we present the results of tapering opti¬ 
mization for the cases of two different superconducting 
technologies, NbTi and Nb3Sn and assess the possibil¬ 
ity of achieving TW levels of power in a 140 m undula¬ 
tor with periodic break sections. We study the impact 
of changing the undulator section length on the perfor¬ 
mance for both NbTi and NbSSn. We then present GEN¬ 
ESIS simulation results for the NbTi case with the opti¬ 
mal choice of section length. Finally we discuss the prob¬ 
lem of particle detrapping as a major source of perfor¬ 
mance degradation in tapered FELs, and propose some 
ideas to improve on current tapering designs by increas¬ 
ing the trapping and the extraction efficiency in the ta¬ 
pered section of the undulator. 

The physical system studied is a 140 m undulator com¬ 
posed of a 20 m SASE section followed by a a 120 m ta¬ 
pered undulator with sections 1 to 3 m in length. The 
SASE and tapered sections are separated by a self seed¬ 
ing chicane which delivers a 5 MW monochromatic seed 
at a photon energy of 4keV. We assume the undulator 
modules to be separated in all cases by 0.5 m break sec¬ 
tions where we install the focusing quarupoles in a FODO 
configuration. The system parameters for both NbTi and 


TABLE I: GENESIS Simulation Parameters 


Parameter Name 

NbTi 

NbSSn 

Electron Beam: 

Beam Energy Eq 

7.2 GeV 

6.8 GeV 

Energy Spread ge 

1.5 MeV 

1.5 MeV 

Peak Current Lpk 

4 kA 

4 kA 

Normalized Emittances ex,nley,n 

0.4/0.4 pL m 

0.4/0.4/i m 

Average /3 function (/3) 

12 m 

12 m 

Undulator: 

Undulator Period Au, 

20 mm 

18 mm 

Undulator Parameter (RMS) a™ 

2.263 

2.263 

Magnetic Gap g 

7.2 mm 

7.2 mm 

Integrated Quad Field Bq 

4.5 T 

4.5 T 

Radiation: 

Photon Energy E-, 

4keV 

4keV 

Peak radiation power input Pseed 

5 MW 

5MW 

Seed laser size Ur 

31 ^ m 

31 /i m 

Rayleigh Range Zr 

10 m 

10 m 

FEL: 

FEL parameter p 

7.67 xlO"'^ 

7.43 xlO"^ 

FEL 3-D gain length Lg~^ 

1.25 m 

1.16 m 

Fresnel Parameter Ed 

8 

8.6 


Nb3Sn are described in table 1 where the beam energy is 
chosen in each case to generate 4keV photons. 

The tapering optimization method used is the one de¬ 
scribed in Ref. m with a transversely parabolic and 
longitudinally uniform electron beam distribution. The 
advantages of using a transversely parabolic beam as op¬ 
posed to a Gaussian in a tapered X-FEL are described 
in detail in Ref. m- The quadrupole gradient is kept 
constant for simplicity and is set to achieve an average 
P function of Pav =12 m throughout the undulator. The 
effect of varying the electron beam size by changing the 
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FIG. 1: Power evolution comparison for different undulator 
section lengths L^ec separated by periodic 0.5 m breaks. The 
simulations are time independent and use the NbSSn param¬ 
eters listed in table 1 starting after the self seeding chicane 
(z=0). The transverse beam size is normalized to the input 
beam size CTe.o and the undulator parameter is normalized to 
the initial value a^o- 




FIG. 2: Time independent GENESIS simulation results for 
the optimized taper profile and evolution of the FEE radia¬ 
tion power for NbTi and NbSSn superconducting undulators. 
We also plot for comparison the case of a permanent maget 
undulator of 2.6 cm period and undulator parameter = 
1.77 as discussed in the text. 


overall extraction efficiency. Thus we determine that to 
obtain a reasonable filling factor and large output power 
a 1.5 m section length is the optimal choice for the un¬ 
dulator design. 


quadrupole focusing strength will further increase the 
output power as discussed in Ref. m however this op¬ 
tion will not be examined in this work. 


II. SIMULATION RESULTS 
A. Section length study 

In this section we investigate the effect of the undulator 
section length on the tapered FEL performance. Since 
the ideal tapering profile is a smoothly varying function 
of z we expect that a longer section length limits the 
performance since shorter sections better approximate a 
continuous magnetic field profile. Furthermore the trans¬ 
verse beam envelope oscillations increase for larger sec¬ 
tion lengths A/32//3^„ = {j3avL)/{Plv - further de¬ 
grading the FEL performance as shown in Eig. 1 for the 
continuous tapering case. We find that for the same av¬ 
erage beta function /3av as the section length becomes 
larger than the gain length the maximum output power 
is significantly reduced as illustrated in Eig. 1. Undula¬ 
tor sections of 1-2 m show a marginal improvement in the 
output power when going from the continuous to the step 
taper which can be attributed to the additional phase 
slippage of the electrons with respect to the ponderomo- 
tive wave in the step-wise case vs the continuous case. 
On the other hand in the 3 m case the maximum power 
decreases from 1.3 TW to 0.5 TW. This is a result of the 
more pronounced oscillations in the beam /3 function and 
the phase mismatch between undulator sections as the 
discontinuity in magnetic field value is larger for longer 
sections. The 1 m undulator sections are thus optimal 
for minimizing oscillations in the beam size and approx¬ 
imating a smooth taper profile however the filling factor 
for such short sections limits the interaction length and 


B. FEL evolution with 1.5 m undulator sections 

For our optimal choice of section length time, indepen¬ 
dent tapering optimizations produce the magnetic field 
profiles displayed in Fig. 2. In both cases the opti¬ 
mal functional form of the tapering law is very close to 
quadratic a^iz) 0^,0 [l ~ c(z — zq)^] with varying ta¬ 
per strengths c and a start location around zq 10 m 
after the self-seeding chicane. The quadratic scaling can 
be understood from the 1-D theory of tapered FELs [T]. 
The change in energy for a resonant electron is deter¬ 
mined by a z dependent poderomotive gradient: 


dir 

dz 


eE{z)ayj{z) 

mc^ 


sin^n 


( 1 ) 


where for simplicity we assume the resonant phase r 
to be constant, E{z) is the magnitude of the growing 
radiation field and 0 ^,( 0 ) is the RMS magnitude of the 
tapered magnetic held. The EEL will radiate at a wave¬ 
length As if the following z dependent resonance condi¬ 
tion is satished: 


7r(^) = ^ (l + a™(2)^) (2) 

If we assume that the radiative process in the ta¬ 
pered section of the undulator is mostly due to coher¬ 
ent emission, with the electric held evolving according to 
E{z) oc z, to satisfy simultaneously Eq. 1-2 the magnetic 
held au,(z) must decrease quadratically with z. With the 
optimal taper prohles for NbTi and Nb3Sn we obtain 
extraction efficiencies of 6.08% and 5.89% respectively, 
and output powers over a factor of 75 larger than the 
saturation power in both cases (see Fig. 2). Applying 
the same tapering optimization for a Permanent Magnet 































































































3 




0.90 

0.85 

0.80 

1 

^ 0.70 
0.65 
0.60 



FIG. 3: Real part of the electron beam refractive index (left) 
and on axis resonant phase (right) for NbTi superconduct¬ 
ing undulator parameters from time independent GENESIS 
simulation. 


FIG. 4: Radiation Phase (on axis) and trapping fraction from 
time independent GENESIS simulations of NbTi. The area 
after 40 m is when the gain is so low that the phase velocity 
of light is almost c thus the radiation phase is constant 


Undulator (PMU) with a 26 mm period, beam energy 
of 6.7 GeV and an RMS undulator parameter au;o=l-77 
we report a ^ 60 % reduction in output power compared 
to the superconducting case with = 1-1 TW at 

the undulator exit. As mentioned previously this is one 
of the motivations for moving from permanent magnets 
to superconducting technology in the design of Terawatt 
level X-FELs such as LCLS-II. We also note that for the 
superconducting case a reduction in the average /3 func¬ 
tion will further increase the efficiency while this cannot 
be done in the permanent magnet case. 


electrons [T]. From Fig. 3 and 4 we see that as the reso¬ 
nant phase increases from zero to ipr = 20° from 2:=10-50 
m, the phase shift in the ponderomotive bucket causes 
~ 20% of the particles to detrap. This can be mitigated 
by introducing phase shifters in the break sections be¬ 
tween 10-50 m and is an effect which will be examined 
quantitatively in future studies. 

D. Time dependent effects. Sideband growth and 
sideband induced detrapping 


C. Limits on optimization due to parasitic effects 

The FEL process is dominated by the effects of re¬ 
fractive guiding and particle detrapping in the tapered 
section of the undulator. In Fig. 3 we show the evolu¬ 
tion of the real part of the refractive index m ( in the 
resonant particle approximation) and the resonant phase 
given by the expressions: 


3fi(n - I) 


, ,2 „2 


UJ 


po ^ &,o [77]/ *-0s 4//^ 

r,l„ I 


Th 2|as 


7fi 


(3) 


sinvl/7,(z)=x^^ (4) 

where y = (2*me*c^/e)(Au,/2As)(l/-\/2[JJ]) is a con¬ 
stant independent of z and the other symbols have their 
usual meaning. From these expressions we see that if 
we want to maintain strong optical guiding, the bunch¬ 
ing must be preserved during the tapered section of the 
undulator. However as the radiation field grows the ef¬ 
fect of optical guiding will decrease as l/|as| inducing a 
self-limiting mechanism on the growth of the field m- 
Furthermore, in order to provide the largest ponderomo¬ 
tive gradient to induce the greatest energy loss in the 
electrons, it is desirable to increase the resonant phase 
throughout the tapered section. Again here we find a self- 
limiting mechanism which limits how much one should 
increase 4 'a, since a larger 41results in a smaller pon¬ 
deromotive bucket creating a tradeoff between the pon¬ 
deromotive potential strength and the number of trapped 


As has been pointed out in previous work [m m nni, 
the performance of tapered hard X-ray FELs is also lim¬ 
ited by particle detrapping due to time dependent effects. 
Shot-noise fluctuations in the electron beam current in¬ 
duce modulations in the radiation beam longitudinal pro¬ 
file. This can cause particle detrapping from the top 
and bottom of the ponderomotive bucket as the modu¬ 
lated radiation field slips through the beam and the par¬ 
ticles experience fluctuations in the bucket height. Fur¬ 
thermore, the resonant interaction between the FEL pri¬ 
mary wave ko and the electron synchrotron motion drives 
a sideband instability which amplifies parasitic modes 
at wavenumbers fco ± Hg where Hg is the synchrotron 
wavenumber: 


(5) 

TTleC^kw 

The growth of the sidebands not only increases the 
bandwidth of the FEL signal but also induces further 
particle detrapping csi- 

Time dependent effects are analysed by simulating the 
beam parameters in table I for a longitudinally uniform 
bunch of 37.5 fs. The simulation results are shown in Fig. 
5 for a longer undulator of 160 m. The simulation was 
done for a longer undulator to emphasize the peak power 
saturation in the final sections of the tapered undulator. 
From Fig. 5 we can see that there is a difference in peak 
output power (averaged over time) compared with the 
time independent results, namely 1.78 TW in the time 
independent case compared to 1.55 TW after 120 m time 
dependent. The drop can be attributed to the growth of 
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FIG. 5: Time dependent GENESIS simulation results for the 
NbTi case parameters in Table 1. The beam is 37.5 fs long 
with a flattop longitudinal profile. 
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FIG. 6: Growth of the FEE synchrotron sidebands at differ¬ 
ent locations in the undulator. The growth of the sidebands 
intensified in the later section of the undulator where the syn¬ 
chrotron frequency is changing less rapidly (see Fig. 5). 


section of the undulator [TS] |50]. In the final section of 
the undulator (z >I00m) the variation in the synchrotron 
frequency is small due to the saturation of the electric 
field. It might be possible to compensate for this effect by 
introducing a cubic or quartic term to the taper profile. 
This option will also be explored in the future. 


III. CONCLUSION 

We have presented a numerical study of tapering op¬ 
timization for a self seeded hard X-ray FEL with super¬ 
conducting undulator parameters. This work has demon¬ 
strated numerically that peak power levels over 1 TW can 
be achieved in a 120 m undulator with break sections and 
an optimized taper profile for two separate superconduct¬ 
ing technologies: NbTi and NbSSn. We demonstrated 
a 60% increase in output power for the superconducting 
technology when compared with a permanent magnet de¬ 
sign. We have also outlined the effects which limit the 
growth of the power in the tapered undulator and have 
found that particle detrapping is the main obstacle to 
achieving larger extraction efficiencies. The two causes 
of particle detrapping we have discussed are due to an 
increase of the resonant phase after the exponential gain 
regime and sideband induced detrapping due to time de¬ 
pendent effects. To mitigate these two effects we propose 
to study applying phase shifters at the location of ini¬ 
tial saturation and introducing a faster term in the taper 
profile to reduce the sideband growth by changing the 
synchrotron frequency faster in the final sections of the 
undulator. 
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